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GEODESIC GROWTH OF RIGHT-ANGLED COXETER GROUPS BASED 

ON TREES 

LAURA CIOBANU & ALEXANDER KOLPAKOV 


Abstract. In this paper we exhibit two infinite families of trees {TnlnAi? and {T^}n>i 7 on n 
vertices, such that and are non-isomorphic, co-spectral, with co-spectral complements, 
and the right-angled Coxeter groups (RACGs) based on and have the same geodesic 
growth with respect to the standard generating set. We then show that the spectrum of a 
tree is not sufficient to determine the geodesic growth of the RACG based on that tree, by 
providing two infinite families of trees and {S'^}n>ii, on n vertices, such that 

and are non-isomorphic, co-spectral, with co-spectral complements, and the right-angled 
Coxeter groups (RACGs) based on and have distinct geodesic growth. 

Asymptotically, as n —^ oo, each set T^, or 5),, i = 1, 2, has the cardinality of the set of 
all trees on n vertices. Our proofs are constructive and use two families of trees previously 
studied by B. McKay and C. Godsil. 
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“And then he had a Clever Idea. He would 
go up very quietly to the Six Pine Trees ...” 
A.A. Milne, Winnie-the-Pooh 


1. Introduction 

The geodesic growth function of a group G with respect to a finite generating set S counts, 
for each positive integer n, the number of geodesics of length n starting at the identity Iq in 
the Cayley graph of G with respect to S. The geodesic growth series is the formal power series 
that takes the values of the geodesic growth function as its coefficients (see Definition I2.4p . 

The groups that we consider in this paper, right-angled Coxeter groups, or RACGs, are 
known to have a regular language of geodesics with respect to the standard generating sets, 
and therefore rational geodesic growth series (see |12] or [3l Theorem 4.8.3] for proofs of these 
facts). Our goal here is to obtain more specihc data concerning the geodesic growth of RACGs. 
Namely, we are interested in extracting information about the geodesic growth series from the 
dehning graph of the group. It is known that non-isomorphic graphs dehne non-isomorphic 
RAGGs [H], and that non-isomorphic RAGGs can have equal geodesic growth [1]. However, 
we are interested in knowing how much the similarities or differences between two defining 
graphs influence the geodesic growth of the corresponding RAGGs. 

Two RAGGs Gi = G{Ti), i = 1,2, with non-isomorphic dehning graphs Fj, may have 
equal standard growth (see Dehnition 12.41 (2)); this can be determined by computing the /- 
polynomials of the graphs Fj [5l Proposition 17.4.2]. However, their geodesic growth exhibits 
more subtle properties [I], and in general it is not known which graph theoretic conditions 
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completely determine the geodesic growth of a RACG. The examples of non-isomorphic graphs 
Tj defining RACGs with equal geodesic growth in [T] are degree-regular and have cycles. In this 
paper we consider the case when the Tj’s are trees. Although these are some of the simplest 
classes of graphs, we already encounter a phenomenon that shows a great difference between 
the behaviour of the standard growth and that of the geodesic growth. The standard growth of 
each RAGG Gi is determined solely by the number of vertices (or edges) in the respective tree 
Tj, while the geodesic growth can be distinct even for two co-spectral trees (which might have 
co-spectral complements as well). Recall that two graphs are co-spectral if the characteristic 
polynomials of their adjacency matrices are the same. 

In this paper we count and compute with the help of an automaton generating the geodesic 
language in a RAGG based on a tree T. This automaton reflects some of the path information 
from T. In general, a good deal of combinatorial information about T can be extracted 
from its spectrum. However, we encounter two rather different behaviours: on one hand 
we produce (infinitely many) pairs of trees Ti and T 2 which are non-isomorphic, co-spectral, 
with co-spectral complements, and the RAGGs G{Ti), i = 1,2, based on them have the same 
geodesic growth; on the other hand we obtain (infinitely many) pairs of trees 5i and S 2 which 
are non-isomorphic and co-spectral, with co-spectral complements, whose respective RAGGs 
G{Si) have distinct geodesic growth. This shows that the spectrum of a tree alone does not 
determine the geodesic growth of the RAGG based on that particular tree. 

The following two theorems, proved in Sections [3] and [3 respectively, are the main results 
of the paper: 

Theorem 1.1. There exist two families of trees Tf = {T^,r 2 ^,... } and Tf = ...} on 

n > 17 vertiees such that for all i >1: 

(1) TI and Tf are not isomorphie, but eo-spectral, with co-spectral complements, and 

(2) the RACGs G(rl) and G{Tf ) have equal geodesie growth series. 


Theorem o thus answers positively Question 1 in [TJ Section 8]. 

Theorem 1.2. There exist two families of trees = {5|, 5^,... } and = {Sf, 51,...} on 
n > 11 vertiees such that for all i > 1: 

(1) Sj and Sf are not isomorphic, but co-spectral, with co-spectral complements, and 

(2) the RACGs G(5f) and G{Sf) have distinct geodesic growth series. 


We note that McKay showed in m that the trees T/ and Tf in Theorem 11.11 are simultane¬ 
ously co-spectral, have co-spectral complements and co-spectral line graphs. On the contrary, 
the trees S} and Sf from Theorem 11.21 might have line graphs with distinct spectra. Moreover, 
by a result of McKay m, in both Theorem 11.11 and Theorem II.21 the cardinality of each family 
Tn and 5^ tends asymptotically to the cardinality of the set of all trees on n vertices: 


( 1 ) 


card 71 , , card 5* 

-—^ 1 and - 

card Tn card 


as n — 00 , f = 1, 2. 


The two theorems above, together with McKay’s results, lead us to the following conjecture: 


Conjecture 1.3. If two trees are simultaneously eo-spectral, have co-spectral complements and 
co-spectral line graphs, then they have the same geodesic growth. 


Section [J] contains details about the computation of several kinds of geodesics, and their 
numbers, for the trees introduced in Section [3l 
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The authors have created a Python code MontjO [H], which performs the computations 
needed for the results in this paper by using either SAGE standard routines for symbolic 
computation of determinants and rational expressions of growth series, or the Berlekamp- 
Massey algorithm for restoring the rational expression for a growth series from a sufficient 
number of its coefficients, as a more efficient approach. Our Python code constructs finite- 
state automata that accept the geodesic languages in RACGs based on triangle-free graphs (in 
this case, trees), and then proceeds to determine the respective growth series. 

2. Definitions and notation 

Let S be a finite set and S* the free monoid on S. We identify S* with the set of words over 
S, that is, finite sequences of elements of S. We use |.| to denote word length. 

Let G = {S) be a group generated by S. For an element g of G, denote by the word 
length of g with respect to S. Given tc € S'*, we denote by w the image of tc in G under the 
natural projection S* — G. 

Definition 2.1. A word w over an alphabet S is geodesic in G = (S) if |r(;| = \w\s- The set of 
geodesics in G with respect to S will be denoted by Geo{G) or Geo{S). 

Let P = r(G, S) be a simple (no loops, no multiple edges) graph with vertex set P(r) = S 
(or simply V) and edge set E{T) (or simply E), where E C V x V. The RACG based on T is 
given by the presentation 

(s € S I = 1 Vs G S, and {ss'f = 1, V{s, s'} € E). 

R is easy to see that for any two involutions s and s' the relation (ss')^ = 1 implies ss' = s's. 
This leads to another possible presentation for RACGs: (s G S' | s^ = 1, ss' = s's, V{s, s'} G E). 

In the present paper we use the same letters for the vertices of P and the corresponding 
generators of the group G(r). 

The star of a vertex u G P in T, denoted by Str(u), or St{v) if the ambient graph is clear in 
the given context, is the set of vertices in P that are adjacent to v. That is, 

Str(u) = {tc G P I {u, rc} G E}. 

We now describe a finite deterministic automaton that recognises geodesics in RACGs (see 
|T0] for definitions of languages and automata). We define such automata in the standard way, 
as quintuples {Q,T,,5,qo,E), where Q is the finite set of states, S the input alphabet, 6 the 
transition function, qq the initial or start state, and E the set of final or accepting states. The 
following dehnition is a simplified version of Proposition 4.1 in [1]. 

Definition 2.2. Let P = (P, E) be a tree. The deterministic hnite state automaton recognising 
the geodesics in G(r) is A = {Q, S, S, {0}, E), where the set of states is Q = {0} U P U E U {p}, 
p is the unique “fail” state, and {0} is the start state. The input alphabet is S' = P, the set of 
accept states T is 0 U P U E, and the transition function 5 : Q x S ^ Q is given by 

(1) 5{cj, s) = (St(s) n u) U {s}, for s ^ a; 

(2) (5((T, s) = p, otherwise. 

Definitions 2.3. Any set L of words over an alphabet S gives rise to a strict growth function 
/l : N ^ N, defined by 

hin) :=|{WgL| |iP|=n}|. 

^this is not the given name of the code, which would be obviously too posh for such a petty thing, but a 
reference name, which is seemingly good for any Python code. 
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Definitions 2.4. Let G be a group generated by S. Then we define the following: 

(1) The geodesic growth function fceo(G) : N —N is given by 

fceoioir) ■■= fGeo(G,S){r) = |{u^ G S* \ |u;| = \W\s = r}|, 
and the geodesic growth series of G equals 

OO 

G(G,S){t) = y^jGeo(G)ir)f. 

r=0 

(2) The spherical (standard) growth function o'(g,S) : N —>■ N is given by 

<7(0,S)(’") = |{5 e G I = r}|, 
and the spherical (standard) growth series of G equals 

OD 

^G{t) ’■= ^(G,S)it) = y~l’7(g'g)(r) 

r=0 

We recall that the f -polynomial of a graph T is the generating function for the number of 
cliques (that is, complete subgraphs) of size i in T: f{t) = /o + fit + f 2 t^ + ■ ■ ■, where fi 
is the number of z-cliques in T. We consider the empty set to be a clique on zero vertices, 
and therefore /o = 1. We remark that the spherical (standard) growth function of a RACG 
is determined by the /-polynomial of its defining graph O Proposition 17.4.2], We want to 
contrast this to the fact that the geodesic growth function is not uniquely determined by the 
/-polynomial, or even by the spectrum of the defining graph, as the following sections show. 

3. RACGS WITH EQUAL GEODESIC GROWTH 

In this section we prove Theorem 11.11 by giving an explicit construction of the families of 
trees z = 1,2. 

Definition 3.1. We define the coalescence r • cj of two rooted trees r and a to be the tree 
which results from merging r and a at their roots. The tree r • a has as vertex set the union 
of the vertex sets of r and a, and as root the identification of the roots of r and a, as shown 

in Fig. [TJ 

Consider the trees Ti and T 2 , both rooted at 0, first described by McKay in m, as shown 
in Fig. [2l In [13] the following fact is proved: 

Theorem 3.2. Let T be a rooted tree with at least two vertices and with root labelled 0. Then 
the trees Fj = T-Tj, z = 1,2, are not isomorphic, but are co-spectral. Also, their complements 
Fj and line graphs L(Fj), L(Fj), L(Fj), z = 1,2, are respectively co-spectral. 

Now let Gi := G{Ti) and G 2 := G{T 2 ) be the RACGs associated to the trees Ti and T 2 
defined above. First note that since the trees Ti and T 2 are isomorphic as graphs, the groups 
Gi and G 2 are isomorphic. However, Ti and T 2 are not isomorphic as rooted trees. 

Let T be a tree with n vertices. Fix a labelling {0,..., zz — 1} of the vertices of r, and suppose 
that 0 represents r’s root. Define Fi = r • Ti and F 2 = r • r 2 to be the trees obtained as the 
coalescence of r with Ti at vertex 0, and let G(Fi) and G(F 2 ) be the RACGs based on Fi and 
F 2 , respectively. Since Fi and F 2 are non-isomorphic, G(Fi) and G(F 2 ) are non-isomorphic, 
as well. 

The following lemma is the key ingredient of Theorem 11.11 
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Figure 1. Trees r and a with marked roots (on the left) and their coalescence 
r • a (on the right) 



Figure 2. McKay’s rooted trees: Ti on the left, and T 2 on the right 


Lemma 3.3. The groups G'(ri) and G'(r 2 ) have the same geodesic growth series. 

In order to simplify the exposition in the proof of Lemma 13.31 we use the notation below. 

Notation 3.4. Let G be a group with generating set T containing the letter 0. Denote the 
set of words in Geo{T) starting with 0 by Geoo{T), the set of words in Geo{T) ending in 0 by 
Geo^{T), and the set of words in Geo{T) ending and starting with 0 by Geog(r). 

Proof, fof Lemm,a \S.,'^) Notice that /Geo(ri)(^) (respectively /Geo(r 2 )(^)) is equal to the number 
of all words of length r in F^ (respectively, F^) minus the number of those words of length 
r in F| that are not geodesics. We denote the number of non-geodesics by /Geo(ri)(^) (and 
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fceo^ir), r espective ly). Since |ri| = |r 2 |, clearly fGeo(ri){r) = fGeo(r 2 )(.r) if and only if 
/Geo(ri)(’") = /Geo(r2)(^)' 

We now show that /Geo(ri)(^) = /Geo(r 2 )(^) foi’ r > 1. Any word in can be written 
as W1U1W2 ■ ■ ■ Un, where tt* G (r \ {0})* and Wi € T*. A non-geodesic word w in F* belongs to 
one of the following sets (or is of the type), depending on its form: 

A : w contains non-geodesics tij G (r \ {0})* or Wj G T^, or both, for some 1 < i, j < n, or 
B : all Wi and Ui are geodesic on their respective alphabets, and there exists I < j < n 
such that OujO is a subword of w with Uj G (Stri(O) n r)*, or 
C : all Wi and Ui are geodesic on their respective alphabets, and w contains a subword of 
the form sOs, where s G (Stri(O) n r)*. 

Notice that the set of non-geodesics is then A\JBVJC, where An(i?UC') = 0 and BnC / 0. 

The number of words in A depends only on the geodesic growth series of r and Ti, and thus 
it will be equal to the number of words of type A in F^. 

The computations in Section S] show that the following identities hold: /Geo(ri)(^) = 
/Geo(T2)(’’)) fGeoo{Ti){r) = fGeoo{Ti){r), fGeo°(Ti){'>') = fGeo°{T2){'^) /GeoO(Ti)(^) = 

/Geog(T 2 )(^) r > 1, as a result of ([2]). This means that there is a length-presenting bijec- 

tion 00 between Geoo(Ti) and Geoo{T 2 ), i.e. for each geodesic w = Ov & Geoo{Ti) there is a ge¬ 
odesic w' = Ov' = (j)o{w) G Geoo{T 2 ), and |t(;| = Analogously, there is a length-preserving 
bijection (fP : Geo^{Ti) —>■ Geo^{T 2 ), and a bijection 0 q : Geoo(Ti) —>■ Geoo(T 2 ). This means 
there is a length-preserving bijection (j) : GeoQ{Ti) UGeo°(Ti) Geoo{T 2 ) U Geo^{T 2 ) between 
those geodesics starting or ending with 0 , for which we have that 0 |Geoo = ^Oi <^lGeoO = and 
^IgcoO = </'o- By the computations in d^]) we have that /(GeooUGeoO)(Ti)(^) = /(GeooUGeoO)(T 2 )(^) 
by the standard formula for the cardinality of the union of two sets, and since /Geo(Ti)(^) = 
/Geo(T 2 )(^)) we also have that /Geo\(GeooUGeoO)(ri)(^) = /Geo\(GeooUGeoO)(T 2 )(^)- Thus, there is a 
length-preserving bijection 0 between Geo \ (Gcoq U Geo^){Ti) and Geo \ (Geoo U Geo^){T 2 ). 

The bijection 0 can be extended to the set Geo(Ti) of all geodesics on Ti by letting 0(tc) = 
'ip{w) for all w G Geo \ (Geoo UGeo'^)(Ti), and then furthermore extended to Geo(Ti) UGeo(r \ 
{0}) by letting 0(tc) = w for all w G Geo{T \ {0}). 

Then 0 provides a bijection between the words of type B in Fj^ and the words of type 
B in F 2 . To see this, associate to each wiUiW 2 ■ ■ - Un the word 0(rci)0(tti)... 4){wn)4>{un) = 
4){wi)ui... (p{wn)un- Then Wi ends in 0, rcj+i starts with 0, and Ui G StT-(O)* if and only if 
(i){wi) ends in 0, 4){wi+i) starts with 0, and 0(tti) = Ui G Str(O)*, by definition. 

It is immediate to see that 0 also provides a bijection between the words of type G, and 
between the words of type i? n G in Fi and F 2 , respectively. Thus, there is a length-preserving 
bijection between the words of type AuBuC {i.e. all non-geodesic words) in Fi and F 2 . This 
concludes the proof of the lemma. □ 

Proof. (of Theorem \1.1\) Let = {ti , T 2 ,... } be the set of non-isomorphic trees on k > 2 
vertices. For f = 1, 2 the two families of trees = {r-Tjlr G T^}, /c = n — 15, satisfy Theorem 

m 

From Theorem 1,1.21 we already know that Fi = r • Ti and F 2 = r • T 2 are co-spectral, for any 
T G Tfc. By Lemmathe groups G(Fi) and G(F 2 ) have the same geodesic growth series. □ 

Theorem 3.5 (Lemma 4.3 in [l3]). Letpi{n), i = 1,2, be the proportion of trees on n vertices 
that have Ti as a limb. Then pi{n) = P 2 {n) for all n and \\m.n^ooPi{n) = 1. 
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From Theorem 13.51 we obtain that —>■ 1, as n —oo. 

4. Computing the numbers of special geodesics 


In this section we prove, by concrete computations, the equalities between the numbers of 
special geodesics required in the proof of Lemma 13.31 

Let Ti and T 2 be as in Figure [2l and recall Notation 13.41 For each group Gi = G{Ti), 
i = 1,2, we construct a finite automaton Ai accepting the geodesic language of Gi, as described 
in Definition 12.21 and with the help of Ai we compute the growth series 

7Gi(t) = fGeoiTi)ir)t\ 07Gi(i) = Y fGeoo{Ti)ir)t\ 
r>0 r>0 

°7Gi(t) = Y /GeoO(Tl)(^■)^^ and o7Gi(t) = Y fGeoO{Ti)ir)t'’. 

r>0 r>0 

We also compute the analogous growth series 7 g2(^)) o 7G2(^)) °7G2(^) and o7g2(^) G 2 , 

where the series coefficients are given by the sequences Geo{T 2 ){r), Geo{T 2 )o{r), Geo{T 2 )^{r) 
and Geo{T 2 )Q{r), respectively. 

Our computations, which we elaborate upon below, show that 
(2) 7Gi(t) = 7G2(i), 07Gi(t) = 07G2(i), °7Gi(i) = °7G2W and o7Gi(i) = o7G2 W- 


These identities prove the equality of the corresponding numbers of geodesics. 


Let A be a deterministic finite-state automaton with accepting states qi, i = 0,1,..., N, 
where qq is the start state and the “fail” state is denoted by q. Let M = M{A) be the 
transition matrix of the automaton A. Computing the generating function 7 a( t) of A is a 
standard technique, and the formula for 7A(i) is 


(3) 


Mw 

^ det(/ - tM) ’ 


where I is the N x N identity matrix, and e = (1,0,..., O)'^ and w = (1,1,..., 1)^ are two 
vectors in Z'^. 


Now let Ai be the deterministic finite-state automata accepting the language of geodesics 
in Gi, i = 1,2, and let Mi = M{Ai) be the transition matrix of Ai. Then the geodesic growth 
series 7Gi(i) is a rational function (determined by the equality ([3])), and the coefficients of its 
numerator and denominator can be easily computed, see [7]. The Python code Monty m 
may perform the above computation either by finding the symbolic determinant det(/ — tMi) 
(usually slow), or by applying the Berlekamp-Massey algorithm (a faster one). This Python 
code starts by creating the finite-state automaton Ai, given a triangle-free graph Ti (in this 
case, a tree), and then proceeds to determining 'yait)- 

Since Gi and G 2 are isomorphic, the equality 7 gi (t) = 7 G 2 (^) is immediate. Here we provide 
an explicit formula for these identical growth series. The output of Monty for both trees Ti 
and T 2 consists of two finite-state automata Ai which are isomorphic, as one only renumbers 
the vertices of T 2 in order to obtain Ti. Each Ai has 32 states and 466 transition arrows, and 
the corresponding growth series are 

7g^(t) = = (1 + t)(i + 2f - 2t^ - 4t^ - t^){l + 5t + lOf^ -h 9t^ - - 26t^ - 34t® 

- 22f + 4t^°)(l -8t- 85f^ - 243f^ - 222f^ 332f® 1194t® -b 1349t^ -b 132t® 

- 1510t® - 2008t^° - 1088f^^ -b 28t^^ + 359t^^ -b 170f^^ -b 15t^^)~^. 
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Now we compute the functions o7Gi(^); ^ = 1)2, which will turn out to be equal, as well. 
However, in this case the equality is not known to hold beforehand: even though the groups G\ 
and G 2 are isomorphic, the image of vertex 0 in Ti under the canonical isomorphism does not 
correspond to vertex 0 in T 2 . We shall explicitly compute the generating function o7Gi(l) for 
the number of geodesic words starting with 0 accepted by each Hj, i = 1,2, and then compare 
these functions. Suppose that the start state of each Ai is qo = 0. Let the corresponding 
transition function 6i be so that 5i{qo,0) = = {0}. If a word w labels a path from to 

an accept state, then the word Ow is a geodesic word in Gi starting with 0. Thus, we have to 
compute the generating function for the number of words starting at g^. and ending 

at an accept state: 

Mi w 

= det{I-tM.y 

where e = (0,..., 0,... , 0)"^, and w = {1,1,... , l)"^. Then we use the fact that o7Gi(0 = 

/Cj 

t • oaciit)- By symmetry, we get °7Gi(i) = olGi{t), i = 1,2. 

By using Monty we obtain 

o7Gi(1) = o7G2(0 = ^(1 + ^)(1 + 2t — 2t^ — — t^){l + At + — 9t^ — 

+ f - 3t® - 3t®)(l -8t- - 243t^ - 2221^^ + 332t® + 119Af + 1349t^ + 132t^ 

- 1510t® - 2008t^° - 10881^^ + 28t^^ + 359t^^ + 170t^^ + 


Finally, it remains to compute the growth series o7Gi(l)- 

There are two kinds of words forming disjoint subsets of the geodesic language of Gi that 
we are interested in: 


(I) the words w starting at the accept state g^. = 5(go,0) = {0} of Ai and coming back 
to it: then the geodesic word Ow starts and ends with a “0” (since there are only 
0-transitions leading to g^. = {0} and no 0-transition coming out of g^J, by Dehnition 

E2](i). 

(II) the words w starting at the state g^. = {0} and ending at a state g^. such that 5(g/., 0) = 
Qrrii 7 ^ Qki and g^^ is an accept state: then the word OtcO will be a geodesic word starting 
and ending with a “0”. 

Let QQ:Gi(t) be the generating function for the words of type (I), and o/3Gi(t) be that for the 
words of type (II). Then, o7Gi(l) = t' oOiGiit) + ' o/^Giit)- We have that 


0«Gi(t) 


with e = w = {0,..., 0,^^, 0,..., 0)^. 

ki 

Analogously, 


]PGi{t) 


Mi w 

det(/ — tMi) ’ 


Mi w 

det(/ — tMi) ’ 


with e = (0,..., 0,^^^, 0,..., 0)^, and w having a “1” at position li for all states g/. as 
described above, and zeroes at all other places. 
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By using Monty we obtain 

= o7G2(i) = t • (1 - 5t - 94^2 - 374t^ - 456t'‘ + 955t® + 4275t® + 5652t^ - 1617t® 

- 16773t® - 24255t^° - 7337t^^ + 26583t^^ + 45100t^^ + 26181t^^ - 12789t^® - 34553t^® 

- 24957t^^ - 3147t^® + 8130t^® + 6288t2° + 1398t2^ - 458t^‘^ - 284:t‘^^ - 24t2^) • (1 + 3t 

-h 2t‘^ - 3t^ - 9t^ - 8t^ + + 3t® - -8t- 85t‘^ - 243t^ - 222t^ + 332t^ + 1194t® 

+ 1349t^ + 132t® - 1510t® - 2008t^° - 1088t^^ + 28t^'^ + 359t^^ + 170t^'‘ + 

Thus, the numbers of special geodesics in Gi and G 2 coincide. 

5. RACGS WITH DIFFERENT GEODESIC GROWTH 

In this section we give a proof of Theorem 11.21 Our construction will be analogous to that 
in Section [3l although we shall use different trees, Si and S 2 , in order to construct the families 
5^ = {r • S'i|r € Tfc}, k = n — 9. Namely, we will use the rooted trees in Fig. [3l 



Figure 3. Godsil’s rooted trees: Si on the left, and S 2 on the right 

We set Fi = r • 5i and r 2 = t ■ S 2 , where r is an arbitrary tree with n vertices labelled 
{0,... , n — 1}. Suppose that 0 represents r’s root. Then the following holds: 

Theorem 5.1 (page 27 of [8]). The trees Fi and r 2 are not isomorphic, but are co-spectral and 
have co-spectral complements. Their line graphs T(Fj), T(Fj) and T(Fj) are not necessarily 
co-spectral. 

Proof. The trees Fi and F 2 are not isomorphic, since Fi has the rooted tree a depicted in 
Fig. |4]as a limb more times than r 2 does. 

The characteristic polynomials of r. Si and Fj {i.e. the characteristic polynomials of the 
adjacency matrices of these graphs) are related by 

(t^^iit) = 4>r{t) ■ (l)Si\{0}{t) + 4>T\{0}{t) • 4>Si{t) - t ■ 0r\{O}(i) • 4>Si\{0}{t), 

according to m Lemma 2.2 (i)]. Given that the trees Si and S2 are isomorphic as graphs 
(though not as rooted trees), we obtain that Fi and r 2 share the same characteristic polynomial 
{i.e. are co-spectral). 

The fact that Fi and r 2 are co-spectral follows from the above and [131 Theorem 3.1 (i)]. 
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Figure 4. The rooted tree a with root 0 


Now, by letting r be the line graph on two vertices we obtain that T(ri), T(ri) and T(ri) 
have different characteristic polynomials from L(r 2 ), LiJ' 2 ) and L(r 2 ), respectively. □ 

Lemma 5.2. Let t be a rooted tree with at least two vertiees, and letTi = t- Si and r 2 = t-S' 2 . 
The RACGs G(ri) and G{T 2 ) have distinct geodesic growth series. 

Proof. We will prove that /Geo(ri)(10) 7 ^ /Geo(r 2 )( 10 ). We use the same argument as in the 
proof of Lemm a [3^ that is, fGeo{ri){r) = /Geo(r 2 )(’’) if and only if /Geo(ri)(’’) = /Geo(r 2 )(’’): 
where /Geo(ri)(^) denotes the numbers of non-geodesics of length r, r > 1. In this proof we use 
the result of our explicit computations with Monty, which shows that /Geoo(ri)(^) = /Geoo(r 2 )(^) 
for r < 8 , but /Geoo(ri)( 8 ) = /Geo 0 (ri)( 8 ) = 8919523 for Gi and /Geoo(r 2 )( 8 ) = /Geo 0 (r 2 )( 8 ) = 
8919522 for G 2 . 

Any word in T* has the form W 1 U 1 W 2 • • .Un, where Uj G (r \ {0})*, Wi G S^, and Wi,Ui 
non-empty except for perhaps wi and Un. A non-geodesic word w in r| belongs to one of the 
following sets (or is of the type), depending on its form: 

(A) it either contains non-geodesics Wi G (t\{ 0})* or Uj G S'*, or both, for some 1 < i, j < n, 
or 

(B) all Wi and ui are geodesic on their respective alphabets, and there exists 1 < j < n 
such that OujO is a subword of w and Uj G (Stri(O) fl r)*, or 

(C) all Wi and Ui are geodesic on their respective alphabets, and w contains a subword of 
the form sOs, where s G (Stri(O) n r)*. 

For the remaining discussion we only consider words w of length 10, and call the number 
of non-empty subwords Ui or Wi the syllable length of w. The number of words of type (A) 
depends only on the geodesic growth series of r and S'!, and thus it will be equal to the number 
of words of type (A) in r 2 . Notice that if the syllable length of tc is < 2 there are no words of 
type (B) or (C), so we obtain the same numbers of words in Ti and r 2 . If the syllable length 
of tc is > 4 then all Wi,Ui are shorter than or equal to 6, and our computations show that the 
numbers of special geodesics of length less than or equal to 6 in Fi and r 2 coincide. Thus, a 
discrepancy may appear only when the syllable length is 3 or 4. The words of type (B) or (C) 
with syllable length 4 have the form w = W 1 U 1 W 2 U 2 (or w = U 1 W 2 U 2 W 3 ), where wi ends in 0 
and W 2 starts with 0, or W 2 = 0. But in this case |r(; 2 | < 7, and our computations show 
that /Geoo(ri)(?’) = /Geoo(r 2 )(^)> up to r = 7. 
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Next, we consider the non-geodesics w of syllable length 3. If they have the form w = uiWiU 2 , 
then wi = 0, and the numbers of geodesics Ui is the same since they are all written over the 
same alphabet determined by the tree r, so no discrepancy in the numbers of non-geodesics 
occurs. Thus, it remains to count the words of the form w = wiUiW2-, where |iui| < 8, wi ends 
in 0, |rti| > 1, € (Stri(O) Hr)*, and W2 starts with 0, \w2\ < 8. Again, if |iUj| < 7 we obtain 

the same numbers of non-geodesics. However, a discrepancy occurs when |iui| = 8 or \w2\ = 8. 

Indeed, the number of such words is ‘2{fGeoo{ri){^) deg.,-(0)) in T* and 2 (/Geoo(r 2 )(^) ‘^6gr(0)) 
in r^. By using Monty, we obtained /Geoo(ri)( 8 ) = /Geo 0 (ri)( 8 ) = 8919523 for Gi, but 
/Geoo(r 2 )(^) ~ /Geo 0 (r 2 )( 8 ) = 8919522 for G 2 , so the numbers of non-geodesics in these two 
groups are distinct. □ 

Now we can finish the proof of our second main result. 

Proof, (of Theorem M.^) This is a straightforward consequence of Theorem 15.II and Lemma [5.21 

□ 
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Appendix 


In this section we give the Python code “Monty” that we used in our computations, with 
comments and remarks. Its on-line copy m can be downloaded as a SAGE worksheet from 


the second author’s web-page. 


We begin by dehning the automaton A that recognizes the language of geodesics of a given 
RAGG G whose defining graph F is given as input. 


def Automaton(t): 

# takes a triangle-free graph, returns the automaton recognising 

# the resp. RACG as a digraph 

CliqueComplex = t.clique_complex(); C = list(); 
for s in CliqueComplex.faces().values(): 

for f in s: C.append(set(f)); 
n = len(C); a = DiGraphO; a.add_vertices(range(n)); 
for i in range(n): 

for V in t.vertices(): 

if not(C[i].issuperset([v])): 

St = set(t.vertex_boundary([v])).union([v]); 
d = set([v]).union(st.intersection(C[i])); 
k = C.index(d); 
a.add_edge((i,k)); 

a.set_vertices({i : C[i] for i in range(n)}) 
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return a; 

Given the automaton A, we then compute the growth function of its accepted language (in 
this case, the geodesic language for G with respect to 5). 

def GrowthFunc(a): 

# takes a geodesic automaton, returns the resp. (geodesic) growth function 

am = a.adjacency_matrix(); 
n = am.nrowsO; 

R = FractionField(PolynomialRing(QQ, ’t’)); 

R.inject_variables0; 
m = diagonal_matrix([1]*n) - t^am; 

M = (l/m.det())*m.adjoint(); 
ee = [0]*(n-1); 
ee.append(1); 
e = vector(ee); 
w = vector([1]*n); 
func = e*M*w; 

return func.numerator().factor()/func.denominator().factor(); 

A clique c in the dehning graph F of the RACG G = G(r) corresponds to a state qc in the 
automaton A. We need the following auxiliary function in order to determine the index of Qc 
represented as a vertex of the digraph A (the automaton) created by the procedure Automaton. 

def Index(a, c): 

# takes a geodesic automaton ‘a’, a clique ‘c’ in the resp. defining graph, 

# returns the vertex of the automaton ‘a’ corresponding to ‘c’ 

1 = None; 

for V in a.verticesO : 

if a.get_vertex(v) == c: 

1 = v; 
return 1; 

Below we compute the growth function oa{t)G as described in the proof of Theorem ll.il 
def GrowthFuncStartO(a): 

# takes a geodesic automaton, returns the growth function for geodesic words 

# starting at the state q, where $\delta(Start, 0) = q$ 

am = a.adjacency_matrix(); 
n = am.nrowsO ; 

R = FractionField(PolynomialRing(QQ, ’tO); 

R.inject_variables(); 
m = diagonal_matrix([1]*n) - t*am; 

M = (l/m.det())*m.adjoint(); 
ind = Index(a, set([0])); 
ee = [0]*n; 
ee [ind] = 1; 

e = vector(ee); 
w = vector([1]*n); 
func = e*M*w; 

return func.numerator().factor()/func.denominator().factor(); 
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Now we define a function that takes as input the geodesic automaton A for a RACG G = 
G(r), a list of cliques I = [cq, ci,..., c^] in the respective defining graph T and returns the 
growth function for geodesic words starting at the state q = 6(go,0) that bring A to any of the 
states described by the cliques in 1 . 

def GrowthFuncStartOEnd(a, 1) : 

# takes an automaton and a list of cliques ‘1’ in the resp. defining graph 

# as input, returns the growth function for geodesic words starting at 

# state q, where $\delta(Start, 0) = q$, and ending at any of the states 

# corresponding to cliques in ‘1’ 

am = a.adjacency_matrix(); 
n = am.nrowsO ; 

R = FractionField(PolynomialRing(QQ, ’t’)); 

R.inject_variables(); 
m = diagonal_matrix([1]*n) - t*am; 

M = (l/m.det())*m.adjoint(); 
ee = [0]*n; 

ind = Index(a, set([0])); 
ee [ind] = 1; 

e = vector(ee); 
ww = [0]*n; 
for c in 1: 

ind = Index(a, set(c)); 
ww[ind] = 1; 
w = vector(ww); 
func = e*M*w; 

return func.numerator().factor()/func.denominator().factor(); 

By using a suitable list of cliques I we can compute the functions Qa(t)G' and o(3{t)Q. Namely, 
in the proof of Theorem ll.il we find 

# the function ${_0~0}\alpha(t)_{G_l}$ 
aOOl = GrowthFuncStartOEnd(al, [[0]]); 

# the function ${_0~0}\alpha(t)_{G_2}$ 
a002 = GrowthFuncStartOEnd(a2, [[0]]); 

# the function ${_0~0}\beta(t)_{G_1}$ 

bOOl = GrowthFuncStartOEndCal, [[1], [1,3], [1,4], [2], [2,5]]); 

# the function ${_0~0}\beta(t)_{G_2}$ 

b002 = GrowthFuncStartOEnd(a2, [[1], [1,3], [2], [2,4], [2,5]]); 

The list 1 = [[1] , [1,3] , [1,4] , [2] , [2,5]] above contains the cliques of Ti corre¬ 
sponding to the accept states q of the geodesic automaton Ai for Gi = G'(ri) such that 
for a state p. The list 1 = [[1], [1,3], [2], [2,4], [2,5]] contains the 
cliques of r 2 with analogous properties, corresponding to the states of the geodesic automaton 
A 2 for the RACG G 2 = G(r 2 ). 

The above described Python procedures are also used to perform the computations in the 
proof of Theorem 11.21 

The on-line version of Monty m contains a variation of the GrowthFunc procedure, called 
GrowthFuncBM, that uses the Berlekamp-Massey algorithm for faster computing. 
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